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, Abstract. The Dirac equation for an electron in the central Coulomb field of a point- 

in ' like nucleus with the charge greater than 137 is considered. This singular problem, to 

which the fall-down onto the centre is inherent, is addressed using a new approach, 
basing on a black-hole concept of the singular centre and capable of producing cut-off- 
free results. To this end the Dirac equation is presented as a generalized eigenvalue 
boundary problem of a self-adjoint operator. The eigenfunctions make complete sets, 
orthogonal with a singular measure, and describe particles, asymptotically free and 
CNJ ■ delta- function- normalizable both at infinity and near the singular centre r = 0. The 

barrier transmission coefficient for these particles responsible for the effects of electron 
absorption and spontaneous electron-positron pair production is found analytically 
as a function of electron energy and charge of the nucleus. The singular threshold 
behaviour of the corresponding amplitudes substitutes for the resonance behaviour, 
O 1 typical of the conventional theory, which appeals to a finite-size nucleus. 

<D ' 

^ ■ PACS numbers: 02.30Hq, 03.65Pm, ll.lOJi 

> 
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1. Introduction 

The radial Dirac equation for an electron in the Coulomb potential 

- aZ/r, (1) 

of a point-like nucleus, where a is the fine-structure constant a = 1/137 and Z is the 
nucleus charge (for formal purposes negative values of Z will be also included into our 
consideration), is a set of two first-order differential equations pQ, [2|. It is convenient 
to write it in the following matrix form 

£#(r) = (e + ^A *(r), re(0,oo), (2) 

„ d k , N 

C = -K7 2 — + -0i + ma 3 . (3) 
dr r 

Here the wave function is the two-component spinor 
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<7j are the Pauli matrices: 

ai= { l J)' iCT2= (°i J)' a *=(l -l)' (5) 
m is the electron mass, and e is its energy. The quantity k is the orbital momentum 

k = -(Z + 1) for J = Z + 

k = / for j = 1- -. (6) 

In what follows we confine ourselves to the lowest orbital state k = — 1. Solutions to 
equation (J2J) are known pQ, [2] in terms of confluent hyper geometric functions. 

We are most interested in the supercritical case ct\Z\ > 1, of which the falling to 
the centre is characteristic. One can exclude one component from (j2J) to reduce it to a 
second-order differential equation (see e.g. [T], [2], 0, [I])- Then the Coulomb term (JTJ) 
in (J2J) gives rise to the inverse-radius- squared singular term —(aZ/r) 2 in the potential, 
which is negative (attraction) irrespective of the sign of aZ. This illustrates the presence 
of the fall-down onto the centre phenomenon inherent to the Dirac equation (J2J) with 
ct\Z\ > 1. This case is of direct physical importance as introducing the effect [E], [Z|, 
[3], jH] of spontaneous electron-positron pair creation by an overcharged nucleus created 
for a short time in heavy ion collisions (see the review [0]), and also the effect of strong 
absorption of electrons by this nucleus - see our previous paper [TUj . (Attribution of 
absorbing quality to a singular centre in other problems and within other approaches 
may be found in [TT], [T2j .) Unlike [TQJ , we find it now more appropriate to be dealing 
directly with the set of two first-order differential equations ©• 

To handle the problem, with the supercritical case included, we define £ (jHJ) as a 
self-adjoint operator, with the Dirac equation (j2J) presenting a generalized eigenvalue 
problem for it, and identify physical states with its eigen- vectors. Stress, that L is 
not the Hamiltonian, defined as Ji = £ — —. There are many ways to present a 
given differential equation as a generalized eigenvalue problem for different operators. 
Our choice (j3J) of C is mainly dictated by the demand that the term, responsible for the 
falling to the centre, be kept in the r.-h. side in eq.(0), to later show itself as a singularity 
of the measure, serving the scalar products in the space, spanned by solutions of the 
generalized eigenvalue problem. In other words, our procedure partially extracts the 
singularity from the interaction and places it into the measure. The singularity of the 
measure offers the physical concept of the singular centre as a sort of a black hole, 
emitting and absorbing particles, which are free near the centre and belong to the 
continuum of eigenvectors of C. 

This is different from what may be obtained using the self-adjoint extension of 
the Hamiltonian within the von Neuman technique, resulting in a discrete spectrum, 
unlimited from below, for the Schrodinger equation case jl], [T^], and in a discrete 
spectrum, condensing near the infinitely-strong-binding point s = —m, for the Dirac 
equation case 
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The present work is an extension to the Dirac equation of analogous treatment, 
developed earlier [H], [T3] for the (second-order) Schrodinger equation with singular 
potential with inverse-square singularity. We find that the extension to the relativistic 
case is straightforward, with the only complication due to the spinorial character of the 
Dirac equation. This is the lack of positive definiteness of the measure in the definition 
of scalar products for negative energy and positive charge. It is shown, however, that 
one can easily restrict oneself to vectors with a definite sign of the norm within one set, 
and also exclude zero-norm vectors by imposing a certain easy-to-observe convention 
concerning the way the ends of the interval (r^, ru), where equation (J2J) is defined - first 
taken finite - tend to their final positions in the singularity points r = and r = oo. 

In Section El we set the generalized self-adjoint eigenvalue problem associated with 
the Dirac equation in a two-side-limited box r^ < r < ru, r-^ — > 0, ru — > oo with 
zero and (anti)periodic boundary conditions, derive the orthogonality relations with a 
singular measure and point a family, labelled by the ratio R = aZ/e, of Hilbert spaces, 
spanned by solutions of the generalized eigenvalue problem, including the ones that are 
5-function-normalizable and correspond to particles, free both near the origin r = 
and near r = oo. In Section El we describe the coordinate transformation £(r) that 
reduces the generalized eigenvalue problem to the standard one by mapping the origin 
r = to infinitely remote point £ = ±oo (depending on the sign of the energy). In 
the new coordinate £ the free character of the wave function behavior near the singular 
point r = becomes explicit, and its normalization to 5-function becomes standard. 
Referring to the one- dimensional barrier problem on the infinite axis £, to which the 
Dirac equation is reduced in our procedure, we find in Section |U the reflection and 
transmission coefficients for the waves, reflected from and transmitted to the singularity 
in r = 0. It is important, that in the Dirac equation, in contrast to the Klein-Gordon 
case, the function, inverse to the transformation £(r), is two- valued in the negative- 
energy continuum. This fact leads to the change of identification of the incident and 
reflected waves near r = as compared to the positive-energy domain, and results 
in absence of the superradiation phenomenon, described in ^Hj, JZj- In Section El 
the transmission coefficient is interpreted within the Dirac-sea picture as probability of 
absorption of electrons by the supercritical nucleus in the positive-energy range, and 
as distribution of spontaneously produced electrons - which outgo to the singular point 
r = - and positrons - which outgo to the singular point r = oo. These coefficients 
are studied in detail and plotted as functions of the nucleus charge and electron energy. 
Their singular threshold behavior near the point a 2 Z 2 = 1 is established. In concluding 
Section we comment on a comparison with the known results jH] about the 
spontaneous pair creation, obtained using the finite nucleus size cut-off, and discuss 
the lines, along which the second-quantized theory of unstable vacuum [IB], [H], [13 
could be applied to the supercritical nucleus, basing on the consideration made in the 
present paper. 
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2. Bliss eigenvalue problem 

There are two singular points of the set of differential equations (j2J, r = oo and r = 0. 
The fundamental solutions behave near the point r = oo as 

e Ti P r r TK ( 7 ) 

where p = \/e 2 — m 2 , £ = 9 ^-. 

It will be convenient for our present purposes, to present the well-known facts about 
this singular point in the following words . 

When |e| < m, the growing solution should be disregarded, whereas the other 
decreases and is thus localized in the finite region of the r-space, far from the singularity 
point r = oo. We refer to this situation by saying: the singularity in the infinitely remote 
point repulses the particle. 

When |e| > m, the both solutions oscillate. They are 5-function-normalizable, 
because are concentrated mostly near the infinitely remote singular point. We refer 
to this situation by saying, that this singularity attracts particles. The particles, 
captured by the infinitely remote singularity are asymptotically free in the sense that 
the fundamental solutions and the spectrum do not depend on Z, k (if it is not fixed), in 
other words, on any terms in (j2J) that do not survive in the limit r — > oo, as compared 
to the constant term e, responsible for the singularity in the infinity. The solutions 
with positive and negative signs in the exponentials in (J2J), the incoming and outgoing 
waves in customary wording, may be thought of as ones, emitted and absorbed by the 
infinitely remote singularity. 

Now we turn to the other singular point of the set of differential equations (J2J), the 
one in the origin r = 0, and shall treat it in exactly the same manner. Two fundamental 
solutions to equation (J2J) behave near the origin r — > like 

r ±i7 , (8) 

where 7 = y/a 2 Z 2 — 1. 

When a\Z\ < 1, the growing solution should be disregarded, whereas the other 
decreases as r — ► and is thus localized far from the singularity point. We refer to this 
situation by saying: the singularity in the origin repulses the particle. 

When a\Z\ > 1 (supercritical charge), the both solutions (JBJ) oscillate, which 
is typical of the problems with singular attractive potential. Within a generalized 
eigenvalue problem for the operator associated with the differential expression £ (JHJ), 
formulated below, they are 5-function-normalizable, because the measure characteristic 
of this problem appears to be singular in the origin, this fact providing the necessary 
divergence of the norm. The particle is thus localized mostly near the origin, attracted 
by the singularity. The particles, attracted by the singular center, are asymptotically 
free in its vicinity in the sense that the fundamental solutions and the spectrum do not 
depend on e, m, in other words, on any terms in (JSJ) that do not survive in the limit 
r — > 0, as compared to the singular potential. The solutions with positive and negative 
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signs in the exponentials (jSJ), may be thought of as waves, emitted and absorbed by the 
singular centre. 
With the ratio 



e 



fixed, equation (J2J) can be written either as 
£V( r ) =e ( 1 + - 1 tf(r 



or as 



C^(r) 



(9) 



(10) 



(11) 



Thus, equation (J2J) becomes a generalized eigenvalue problem of the type, where the 
eigenvalue (this is either e for (|10j) or aZ for ([11))) is multiplied by a function of the 
variable r, and not by unity In the case of a second-order equation an analogous 
eigenvalue problem was first studied by E. Kamke [20]- we called it Kamke eigenvalue 
problem when using it as applied to the singular Schrodinger equation in [T3], [TU], - 
but for the set like (j2J) it was considered, according to [20J, much earlier by G. A. Bliss 



The equation, in which the sign in front of the derivative term is reversed 
and, besides, all the matrices are transposed, is called an adjoint equation. 
Equation (j2j) is self — adjoint in the sense that solutions \l/(r) of its adjoint equation 
are linearly expressed in terms of its solutions with the help of a nondegenerate matrix, 
namely: \I/(r) = cr 2 ^(r), det o<i ^ 0. The eigenvalue problem ()10|) or (fTTj) is self— 
adjoint, provided appropriate boundary conditions are imposed. As such, it is sufficient 
to use the zero and periodic (or antiperiodic) boundary conditions at the ends of the 
interval. 

We shall consider equation (fTUj). (fTTj) in four intervals: (0,oo), (r L ,oo), (0,ru), 
( r L, r u)> r L; r u > depending on the regions - called sectors - within which the 
parameters e and aZ may lie, and pass to the limit r^ — > 0, ru — > oo afterwards. 
So, in the end, all the intervals are one interval. The boundary conditions are: 



G(0) = G(oo) = 0, when \aZ\ < 1, |e| < m; 
the interval is < r < oo 



Sector I 



G(0) = G(r v ) = 0, when 
the interval is < r < r\j 



\aZ\ < 1, lei > m; Sector II 



G{ri) = G(oo) = 0, when \aZ\ > 1, \e\ < m; 
the interval is r^ < r < oo 



Sector III 



G(r L ) = ±G(r v ), F(r L ) = ±F(r v ), when \aZ\ > 1, |e| > m; 
the interval is r^ < r < ru Sector IY. 



(12) 



(13) 



(14) 



(15) 



Black-hole concept of a supercritical nucleus 



6 



As long as tl 7^ for a\Z\ > 1 and/or ru ^ oo for |e| > m, there are infinite 
countable manifolds of eigenvalues in sectors II, III, IV, which condense in the limit 
rL — > 0, r\j — > oo to make continua. 

The function e + — does not change sign throughout the interval (0, oo), when 
R > 0, i.e. in the positive energy domain (e > 0) for positive charge Z > 0, and 
in the negative energy domain (e < 0) for negative charge Z < 0. In these regions 
the self-adjoint boundary problem is positively definite. In this case the eigenvalues are 
real and the corresponding solutions make, for each value of R, a complete set of states, 
mutually orthogonal with the measure 



Below we shall consider also negative R and see, what happens in that case. Note, that 
the e + e~- pairs production just occurs in the region R, of sector IY. 

The domains, pointed in (fT2|) . (JT3j) . (fTljl. (fT5]). correspond to what was called sectors 
I, II, III, IV in jXHJ, [2]. Sector II (a\Z\ < 1, \e\> m) corresponds to particles, free at 
infinity in the limit r\j — > oo and repulsed from the centre. Only one of two fundamental 
solutions belongs to L^(0, ru), the space of functions square-integrable with the measure 
(|TBJ) in the interval (0,ru). Sector III (a\Z\ > 1, |e| < m) contains particles, free in 
the origin in the limit — > and repulsed from the infinitely remote point. Only one 
fundamental solution belongs to L 2 Jri^ oo)). In sector IV (a\Z\ > 1, |e| > m) particles 
are free both near the origin and near the infinity in the limit j~l — > 0, ru — > oo. The both 
fundamental solutions are L^(r L ,ru). The (anti)periodic boundary conditions (fT5|) are 
imposed in agreement with the inelastic character of the scattering in sector IV, where 
there is a nonzero current inflow through the outer border r = ru, equal to the outflow 
through the inner border r = rL, and vice versa. As distinct from this, in sectors II 
and III the current is zero, the scattering of particles emitted by the infinity (sector 
II) and by the centre (sector III) is elastic: everything what is emitted by the centre 
(infinity) is reflected back to the centre (infinity). As for sector I (a\Z\ < 1, |e| < m), 
where the infinite countable manifold of hydrogen-like bound states lies in the quadrants 
(0 < aZ < 1, < e < m) and (— 1 < aZ < 0, — m < e < 0), neither fundamental 
solution generally belongs to L^(0, oo). 

To be more explicit, let us - following the standard way - left- multiply the equation 
(110)1 (or (fTTjl ) by the spinor ^*(r), which is the solution of the complex-conjugate 
equation, but with e\ taken instead of e (or aZ\ instead of aZ) and the same R. (Stress 
that all the coefficients and matrices in expression (JHJ) are real and do not depend 
either on e or on aZ). Let us, next, subtract the same product, with \l/(r) and \E^(r) 
interchanged and integrate the difference over r within the limits r\ < r < r<i- In 
agreement with (fT2|) . (JT3j) . (JHJ), (JUJ), in sectors I and II the lower limit should be 
chosen as t\ = 0, whereas in sectors III and IV it is t\ = r^. The upper limit r2 
coincides with ru in sectors II and IV, and is infinite r2 = oo in sectors I and III. Then 




(16) 
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one has 



- i *I(r)«T 2 *(r)|^ = (el - e) jT *J(r)*(r) (l + dr = 

= a(ZJ - Z) jT *J(r)*(r) (± + ^ dr. (17) 

When written explicitly in components this equation looks like: 
(F?(r)G(r)-F(r)Gt(r))\:i = 

= (e\ - e) jf 2 (G(r)G*(r) + F(r)F*(r)) (\ + |) dr = 
= - Z) f\G(r)G\(r) + F(r)F*(r)) f i + -) dr. (18) 



The boundary conditions (|12|). (|14j) . (|13|) . (|15|) provide the vanishing of the left- 
hand side of (fTTj) or (JTSJ) and make the differential expression C a self-adjoint operator. 
By omitting the index 1 in this relation we obtain that for the eigenvalues e and Z of 
the self-adjoint operator C to be real 

e = e\ Z = Z* (19) 

it is sufficient that the function (-^ + -) should have a definite sign, since in this case 
the integral in (fTTj) is nonzero. Such situation occurs when R > 0. 

Otherwise, when R < 0, the zero-norm vectors, satisfying the relation 

r2 ( R\ 

#*(r)#(r) ^1 + —J dr = (20) 

might exist and then the eigenvalues corresponding to it be not necessarily real. 
Also negative-norm solutions might appear with real eigenvalues, and the manifold of 
eigenfunctions of the self-adjoint operator (j3J) might not be complete, according to [20J, 
[21] • In such case it had to be completed to the Hilbert space, as explained in (2*2 . As 
a matter of fact, this scenario can be avoided. 

It follows from IflKj) . (fT2"|) . (fT3|) . (JHJ), (fT5|). that the eigenfunctions with nonzero 
norm obey the orthogonality relations 

n(r)V n (r) 1 + - dr = JV<W (21) 



r , 

Here the integers fc, n label the eigenfunctions belonging to different countable 
eigenvalues e or Z, and the wave functions \l/ relate to a fixed value of R. The norm 
N in every sector, except sector I, gets predominant contribution from the integration 
near the points r = tl and r = r\j. 

Bearing in mind the asymptotic behaviours (J7|) and ©, and the boundary 
conditions ()15|) . we get for the eigenfunctions in sector IV in the asymptotic regime 
of very large upper size of the interval r\j ^> |i?|,m _1 and very small lower size 
r L <C \R\, m _1 - up to a finite factor - that 

N IV x J T J + ^ dr x ru + i? In X (22) 
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Here r is an arbitrary positive dimensional constant, <C r <C r\j. We neglected tl 
as compared to | R In ^ | , and | i? In 1 as compared to r\j . A change of ro does not 
violate eq. (|2"2"j) within its accuracy. 

For positive R the norm ()22|) is positive. For negative R the norm ()22|) /ias a definite 
sign for a definite value of i?, provided either of the inequalities 

-j — r > In — or - — r < m — (23) 

is observed during the limiting process — > 0, ru — > oo. By imposing this additional 
requirement we avoid the appearance of vectors with different signs of the norm within 
the same set of eigenfunctions. The zero-norm vectors of state do not appear either. 
Thus, the problem of indefinite metric is eliminated in sector IV. In other sectors it does 
not appear. 

In sector II, where the eigenfunctions are L^(0,ru), the integration in (J2"T|) near 
r = is convergent, while the upper limit gives a predominant contribution to (J21|) for 
ru ^> \R\,m~ l due to the asymptotic behaviour (|7J) 

iV" x J U ^1 + dr x ru. (24) 

This has a definite sign. In sector III, where the eigenfunction is L^(rL,oo), the 
integration in (}2"Tj) near the upper limit converges, but near the lower limit gives a 
predominant contribution for tl <C |-R|,m _1 due to the asymptotic behaviour (jHJ) 

JV /w x / + (25) 

Jr L V r J r L 

This has a definite sign for a given sign of R. As for sector I, the eigenvalue problem 
(fTU|) . (fTTj) with the boundary conditions (fT2*|) does not have any solutions for R < 0: the 
hydrogen-like bound states only lie in the segments with R > in this sector. In the 
latter region, for the known discrete values of the energy e, depending on the charge 
aZ, which provide the fulfillment of the boundary condition (|T2*j) . the solution decreases 
as exp{— r^Jm? — e 2 } at r — >• oo and as r^ 1- " 222 at r — >• 0. Correspondingly, the norm 

N 1 x J #*(r)tf(r) ^1 + dr < oo (26) 

converges in spite of the singularity in the measure. 

It may be seen from the results of 0, where the Klein- Gordon equation was 
considered, that the corresponding measure would in that case be d/i(r) = + ^) 2 dr, 
which is positively definite, and hence the problem of indefinite metrics does not arise. 

Manifolds of eigenvalues of the problem ()1U|) (or ([11))) may belong simultaneously to 
two or three sectors: to sectors I, II, IV, if \R\ < m _1 , to sectors I, III, IV, if \R\ > m _1 , 
and to sectors I, IV, if \R\ = m~ l . Correspondingly, the resolution of the unity includes 
integration over continua of eigenvectors in all the sectors II, III, IV crossed by the 
straight line R = const in the plane (e, aZ), the states of both signs of the charge and 
energy being simultaneously involved, as well as summation over bound states in sector 
I for R > 0. To conform the boundary conditions (fTB*j) in sector II and (fTHj) in sector 
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IV, what is needed when > R > — m" 1 , we have to choose the left inequality in (J23)) . 
so that the signs of the norms be the same - positive - within one set of eigenvalues. 
Analogously, when R < — m _1 , we have to choose the right inequality in it, so that 
the boundary conditions (jl4j) in sector III and ()15|) in sector IV might be in mutual 
agreement and the sign of the norm be common - negative in this case. At last, in the 
special case R = —rrT 1 the problem is indifferent to the choice of the inequality sign in 
()23|) . The unrestricted growth of the norm ()22j) in the limit r^ — > 0, r\j — > oo gives the 
possibility to replace the r.-h. side of (|2"Tj) by 5-function of (e\ — s) or of [Z\ — Z) . 

The integrals (jUJ), (J21|), play the role of effective sizes of the quantization boxes 
in the corresponding sectors. Once the spacings between the levels prior to the limiting 
transition r^ — > and/or 7\j — > oo are inversely proportional to these effective sizes (cf 
[Hj), we conclude that the densities of states are essentially different in different sectors, 
which is important for physical processes occurring when a change of parameters causes 
transitions between sectors, as in heavy ion collisions. 

Besides the family of Hilbert spaces labelled by the ratio R ©, we can point two 
more families. The standard one is associated with the representation of the Dirac 
equation (J2J) in the form of the eigenvalue problem for the Hamiltonian 

H = C - — , (27) 
r 

the operator, which leaves aside the term of the Dirac equation, singular near r = oo, 
but includes all the terms singular in the origin 

H^(r) = £*(r). (28) 

This eigenvalue problem is well-defined within our procedure, which is standard in this 
case, when there is no singular attractiveness in the origin, i.e. for a\Z\ < 1. The 
families of eigen-vectors, orthogonal for different values of e with the measure dr, are 
labelled by the charge —137 < Z < 137 and consist for each charge of the wave functions 
corresponding to discrete bound states in sector I, — m < e < m, and two continua of 
positive, £ > m, and negative, e < —m, energies in sector II. 

The other one is associated with the representation of the Dirac equation (0) in the 
form of the generalized eigenvalue problem for the operator C — e, which leaves aside 
the term of the Dirac equation, singular near r = 0, but includes all the terms 
singular in the infinitely remote point, 

(C-e)*(r) = — tt(r). (29) 
r 

This eigenvalue problem is well-defined within our procedure, when there is no singular 
attractiveness in infinitely remote point, i.e. for |e| < m. The families of eigen-vectors, 
orthogonal for different values of Z with the measure ^dr, are labelled by the energy 
— m < e < m and consist for each energy of the wave functions corresponding to discrete 
bound states in sector I, —137 < Z < 137, and two continua of positive, Z > 137, and 
negative, Z < —137, charges in sector III. 
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3. Dilated coordinate space 

The generalized eigenvalue problem (jlO)) (or (jllj) ) can be formally reduced to the form 
of the standard eigenvalue problem like 



with the help of the dilatation of the variable 

£(r)= / ( 1 + - j dr = i2In — + r - r . (31) 

The positive dimensional constant ro here is arbitrary. 

If R > 0, the transformation (}3*Tj) is a monotonous function of r, and the function 
r(£), inverse to £(r), is single- valued. In this case the spinor in (J30|) should be understood 
as ^(0 = \l/(r(£)), and the interval, where this equation is defined, is £ G (—00,00). 
The denominator in the second term of (J30|) is nonzero. 

If, however, R < 0, the denominator turns to zero in the point r = —R, and also 
the function r(£), inverse to £(?"), is two- valued. Denote the two branches of the inverse 
function as r a (£) and rb(£): 

< r a (£) < -R, Um < f < 00 

-i? < r b (£) < 00, £ min < f < 00. (32) 

Due to (J31j) the extremum condition (d£/dr) = is satisfied in the point r = —R, hence 
£min — £( — R)- The value of £ min contains arbitrariness due to the arbitrariness of r . 
Its specific choice is not important. Now (fHUJ) becomes two equations 

. d§ a , b (0 r a , b (0 / \ ~ ~ 

d£ r a , b (£) + R V r a,b(0 / 

£min < £ < OO (33) 

for the two spinors ^ a ,b(0 = ^( r a,b(0)> defined in the same interval of £, but with two 
different functions r a (£) and r b (£) that coincide in the end-point r ai b(£ m in) = — -R- Each 
of these equations has a singularity at the end of the interval £ = £ min . The two spinors 
are connected by the continuity relation 

* a (-fl) = * b (-fl). (34) 
The current 

j(r) = M*(r)a 2 ^(r) = F(r)G*(r) - F*(r)G(r) (35) 

is conserved, (dj/dr) = 0, on the interval < r < 00 due to eq. (fTU|) . Under the 
transformation (|3*T|) it turns into two currents j a ,b(0 = j( r a,b(0) = i^a b (0 a 2^a,b(0 
that both conserve, (dj aj b/d£) = 0, as a consequence of eqs. (|3^j) . in the interval 
£min < £ < 00. Eq.(jSU) guarantees the coincidence of the two currents in the singular 
end-point £ = £ m i n . We conclude that the current does not undergo a discontinuity 
when passing the singularity point in equations (}3*3"j) . once the latter are obtained from 
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The transformation (|31J1 maps the point r = oo to £ = oo, and the point r = to 
either positive or negative infinity, depending on the sign of the ratio R (JSJ). In terms of 
the delated variable £ the asymptotic behaviour (JHJ) - up to the coordinate-independent 
factor r^ 17 - becomes 

exp(±i^), |£|-oo. (36) 

This looks like an ordinary free wave in the £-space with the psewc/omomentum (j/R) = 
±V 1 e 1 — R~ 2 . (For the special value R = — m , noted above, the psedomomentum 
coincides with the momentum.) Under the transformation (|31|) the measure ()16|) turns 
into d/i(r) = d£, and the orthogonality relation (J2TJ) now is 

tt£(0*n(0<£ = (37) 

'£(n) 

The norm (|22|) in sector IV becomes 



iV JV = / d£ = £ L + £u, if J R>0 



and 



iV "' ;s/ ^ 1 + ^)) dr ' ({) + /I( 1 + 4 ldrb « ) 



I'L 



.in /"£u 

d£+ / d£x£ L + £u, if i?<0 (38) 

where £u = £( r u) x r u — > 00 an d £l = ~£( r L) x -Rhi( r o/ r L) —* sgnR ■ oo. Thus, in 
the representation of the dilated variable £ it is explicit that eigenfunctions in sector IV 
are normalizable to S- function, when £l+£u tends to positive infinity. This occurs, when 
R > and when — mr x < R < 0. In the latter case £l < 0, but £l + £u —* +oo, once 
the sign > is chosen in (J23J) . as explained in Section 2. On the contrary, for R < — m _1 , 
we agreed above to choose the < sign in (J2HJ), so that £l + £u — * — oo, and solutions 
are normalized to — S(e — e-\). Then the total probabilities should be defined as norms 
taken with the minus sign. If, at last, R = — m _1 , either convention about the sign in 
and, correspondingly, about the total probability does fit. 
Analogously, the norm (|2*5j) in sector III is 



N UI x / d£ = £ L -> sgni? • oo. (39) 

This is normalizable to 5-function for R positive, and to minus 5-function for R negative. 

The contents of this section perfects the demonstration of complete symmetry 
between the inner and outer coordinate spaces, r — > (|£| — > oo) and r — > oo (£ — > oo). 



4. Transmission and reflection coefficients 



In sector IV, free particles emitted by infinitely remote point (waves, incoming from 
r = oo) are partially reflected backwards, but partially penetrate to the vicinity of the 
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origin to become free near it (waves, outgoing to r = 0), in other words - absorbed by the 
centre. Also the inverse process may proceed: free particles emitted by the centre (waves, 
incoming from the origin r = 0) are partially reflected back to the centre, but partially 
escape to infinity (waves, outgoing to r = oo and " absorbed by the remote point" ) . Thus, 
the singular differential equation is no longer a one-particle problem, but corresponds to 
an inelastic, two-channel problem, characterized by a 2 x 2 scattering matrix, the same 
as in (J2|), ([El)- Below we shall find the off-diagonal S'-matrix element squared - the 
transmission coefficient that is responsible for absorption of electrons by the supercritical 
nucleus, provided that the electron energy is above m, and for spontaneous elecrton- 
positron pair production, provided that the energy is below —m. The latter conclusion 
holds true within the concept of the Dirac sea - usually appealed to, when this effect is 
considered in the customary context (see the review jHJ, and the monographs jH]). The 
negative-energy free electron, which at r = oo belongs to the filled Dirac sea, after it is 
scattered off the point-like nucleus, partially transmits to the centre to become free near 
it and leaves an empty vacancy in its stead, interpreted as a positron, free at r = oo. 
Unlike this scenario, the one of refs. [3], [S] specifies that the electron near the nucleus 
is not free, but occupies the deepest bound state. 

In the present section we deal only with positive nucleus charge Z and are in sector 

IV 

aZ > 1, \e\ > m. (40) 
Consider the fundamental solution to equation (j2J), which behaves as a single exponent 



near r = 



M/(r)U ~(2pr) i7 



^ + 0(^-1)^ 



(41) 



where 



p = y/e 2 - m 2 , 7 = Va 2 Z 2 - 1, C = (aZe)/p, |C| > 7 (42) 
are real quantities. Besides, \(\ > 7. The same solution behaves at r = 00 as 

+e^(2pr)<(^ ^1^1)1 V ( 43 ) 

The coordinate-independent coefficients 8 and Q can be determined using the analytical 

continuation yl — a 2 Z 2 = \\Ja 2 Z 2 — 1 of the known exact solution (see PP, [2]) into 
supercritical region a 2 Z 2 > 1. These are 

r(2i 7 + l)e-g«+7) r(2i 7 + l)e-*«-r) 

i(7-C)r(i( 7 -C))' (i-ifC)r(i( 7 + C))' 1 } 
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Calculating the values of the conserved current (|35|) near r = and near r = oo, and 
equalizing the two values with one another we obtain the unitarity relation in the form 
of the identity 



\Q\ 2 = \S\ 2 + p-, (45) 
C -7 

easy to verify explicitly by the substitution of (44). 

In the positive energy domain, e > m (R > 0), in accord with (|36p. the exponential 
in (|4~T|) r 17 ~ exp (i£^) oscillates with the same sign of frequency as the second term in 
(|43|) . Referring to the second term in (|43|) as to the wave, incoming from infinity, and 
to (j4*Tj) as the wave transmitted to the centre, we have to normilize the incoming wave 
to unity. This reduces to division of eq.([45|) over \Q\ 2 to give it the form 

1 = R + + T+, (46) 

with M + and T + being the reflection and transmission coefficients, respectively: 

R+ = W T+ = (c^WF' £>m ' (47) 

As ( > 7, when e > m , one has T + > 0, hence M + < 1, T + < 1. Finally the following 
analytical expression can be deduced for the transmission coefficient ()47)1 from (jUj) 

T + = 1 - jgj! = 1 - e^'« - 1\ = 2 f 1* 2 7, £ > m .(48) 
+ \G\ 2 sinh7r(C + 7) e 27r «+ 7 ) - 1 v ; 

For sufficiently large 7 this becomes 

1 



T 1 



1 _ e -27r( 7 +C) • 

Owing to the two-fold character of the transformation inverse to (|31|) in the negative 
energy domain, the identification of the incoming and reflected waves is there different. 
Now that e < —m (R < 0), in accord with (JHHJ), the exponential in (JIT]) r 17 ~ exp (i£^) 
oscillates with the same sign of frequency as the first term in (J43|) . it is the first term 
that should be referred to as the wave, incoming from infinity, and the second one as 
the reflected wave. Now the normalization of the incoming wave to unity implies the 
division of eq. ()45|) over \S\ 2 . Then eq.([45|) acquires again the form (J4T)j) 

1=M_+T_, (49) 

but this time the reflection and transmission coefficients are 

R - = jJ T - = (^w £< " m - (50) 

In this domain ( < —7, and one has T_ > 0, hence M_ < 1, T_ < 1. Finally 
the following analytical expression can be found for the transmission coefficient in the 
negative energy domain 

T = 1 _ K. = T + = ! _ e 2, T gjjMC±li 
\£\ 2 T+-1 sinh7r(C-7) 

1 - e" 4w7 , , 
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For sufficiently large 7 this is simplified to 

T _ = e 2^(C+7)_ (52) 

The latter expression is valid already for aZ — 1 ^> 1/327T 2 . 

In our preliminary publication JU] we did not take into account the fact that 
the two-valuedness of the coordinate transformation, inverse to pip, affects the 
identification of incident and reflected waves, and took (|4~Tj) as universal expressions 
for the reflection and transmission coefficients valid both in the positive- and negative- 
energy domains. The reflection coefficient R+ ()47jl . when extended to e < —m becomes 
greater than unity, while the transmission coefficient T + (|47j) becomes negative. Facing 
such situation, known as superradiation, we concluded that the spontaneous pair 
production is impossible due to the Pauli ban: electrons, taken from the Dirac sea, 
would have been increased in number after reflected off the nucleus. This is, however, 
forbidden, since all the vacancies in the Dirac sea are filled. We now renounce this 
point of view: the reflection and transmission coefficients in the negative-energy domain 
defined by (JHUj) and (jHTj) obey the inequalities < IR_ < 1, < T_ < 1, the same 
as (|47p. (}4*H|) in the positive-energy domain. In other words, the superradiation does 
not occur in our case of a spinor particle, where it might forbid the spontaneous pair 
creation. 

The following comparison with the known results on scattering of Bose- and 
Fermi-particles off a black hole is in order. According to ^Hj, if the black hole is 
rotating and correspondingly described by the Kerr metrics, the dilation transformation, 
analogous to (}30|) that reduces the corresponding propagation equation to a standard 
form of a Schrodinger-like equation may, in a certain kinematical domain, be two-fold- 
reversible both for electromagnetic and gravitational waves, and the Dirac spin-| field. 
Simultaneously a singularity comes out from beyond the horizon into the coordinate 
region where the differential equation is defined. The general fact is that the current 
(or the Wronsky determinant) changes its sign when crossing this singularity for Bose- 
fields, but does conserve for Fermi-fields. (We observed the same property of current 
conservation across the singular point r = —R in the differential equation f)30|) in Section 
El Also this property holds true for the singular barrier problem associated with the 
second-order Schrodinger-like equation, to which the Dirac equation in the Coulomb 
potential was reduced in j3] and - after a coordinate-dilation transformation - in 
Stress, that the singularity we are referring to is not the original singularity of the 
initial equation, but the one acquired in the course of its transformations.) Therefore, 
the superradiation takes place where there is the current discontinuity in the singularity 
point, i.e. for Bose-particles. 

5. Absorption of electrons and production of electron-positron pairs 

We now turn to more thorough consideration of expressions (|4~H|) and (|5ip . derived above, 
and of their physical implementation. 
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The transmission coefficient (J48|) . when multiplied by an electron flux density, is 
the absorption probability of electrons, incident on the point-like nucleus with Z > 137. 
Substituting ()42j) in (|4"H|) we get the transmission/absorption coefficient as a function of 
e and Z, T + (e, Z). It has the following asymptotic values 

T + (oo, Z) = l- e- 2 ^ 7 < 1, 

sinn 7r(o;z + 'yyZ)) 

T+(m,Z) = 1 - e - 4 ^ 7(z) < 1, 

T + (oo,Z) > T+(m,Z) (53) 

that depend on the nucleus charge Z (see Figure 1). Near the threshold of absorption 
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Figure 1. Transmission/absorption coefficient T+(e, Z) (|48|) . (|42|l plotted in 
logarithmic scale against electron kinetic energy > 1 for three values of nucleus 
charge, from bottom to top, 7 = 0.498, 0.500, 0.502 (aZ = 1.117, 1.118, 1.119). The 
transmissiom/absorption coefficient is close to unity 



7 = [olZ = 1) at the border of sector IV these asymptotic values behave as 

T+(oo, Z)\ aZ ^ ~ 2n v / (aZ) 2 - 1 (1 + cothTr) = 4:.008n^(aZ) 2 - 1, 
T + (m,Z)\ aZ ^ ~ i7ry/(aZy - 1. (54) 

This means that near the very threshold the absorption is very low. On the contrary, 
already for aZ ~ 1.12 (7 ~ 0.5) the value of T + (e, Z) exceeds 0.998 in the whole energy 
range m < e < 00 (see Figures 1, 2). 

The threshold behaviour of (f4"Kj) for any energy is 

T+(e,7)U->i - 2n^/(aZy - 1 (1 + coth^L=), £ > m. (55) 

V e — 7n 

It is also of interest to consider the transmission/absorption coefficient as a function 
of £ and R, since the Hilbert space is formed by solutions with fixed ratio (JHj), according 
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Figure 2. Transmission/absorption coefficient T+(e, Z) l|48|) . I|42|) plotted against the 
nucleus charge for the electron kinetic energy value e — m = 11 • I0~ e m. The curves 
for other energy values are indistinguishable in the scale of the figure from the plotted 
one: the difference with the transmissiom/absorption coefficient for e — m — 100 m 
makes about 10 . 



to Section^ Call T + (e, R) the function, obtained from P5|l by the substitution (R > 0, 
once e > m, aZ > 1) 

e 2 R 



7 = Vem* - 1, C= r~, t (56) 



?71 



Now the asymptotic value 

f + (oo,i?) = l (57) 

is universal for every R. A family of curves T + (e, R) is shown in Figure |3J 

If R > m -1 , the border of sector IV is at e = m, 7 > 0. The values of T + (e, R) at 
this border are 

f+(m, i2) = 1 - e-Wm^^-i^ R > m -i ( 5g ) 

and correspond to crossings of the axis e — m by the family. (Note, that the ordinate 
axis in Figure 3 corresponds to e — m = 10 _5 m _1 , and not to e — m = 0.) 

If R < m -1 , the border of sector IV is at e > m, 7 = 0. The curves in Figure H3 
cross the abscissa axis in the points e = e t hi = R 1 , found from the equation 7 = 0. 
The threshold behaviour of the transmission/absorption coefficient near these points is 



T+(e, R) ~ 2n(2Ry/ 2 V^R—i(l + coth = ), 

e-»e t hr=-R" 1 V 1 — m z R 2 

e > R- 1 > m. (59) 

The absorption takes place, for R fixed, for electron energies exceeding the threshold 
values e t hr = -R -1 - 
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Figure 3. Transmission/absorbtion coefficient T+(e, R) l|48|) . (|56|l plotted in 
logarithmic scale against electron kinetic energy > f for five values of the ratio 
©, from left to right, Rm = 1.003, 1.001, 1.000, 0.999, 0.990,. The curves with 
Rm < 1 (drawn solid) show the threshold behaviour (|59"|l 



In the negative-energy domain e < — m of sector IV, aZ > 1, R < the 
transmission coefficient (JoTj) . when multiplied by the (degenerate) Fermi-distribution 
of electrons in the Dirac sea, which is unity, becomes the distribution of positrons, 
spontaneously produced from the vacuum. These distributions as functions of energy, 
T_(e, R), are presented for the fixed negative ratio R in Figure 4, showing the 
transmission coefficient (|HT|) . with substituted for 7 and ( in it. The asymptotic 
value for large \e\ — » 00 

f _(e, R) ~ e 27r «+^ ~ e -^ Rm+ ^ (60) 

is equal to unity for every R: 

f_(-oo,i2) = 1, ,R<0. (61) 

This means that the total probability of creating a positron with its energy less than 
infinity is 1. 

Within the Dirac-sea picture, we are sticking to, the mechanism, leading to the hole 
distributions like the ones drawn in Figure 4, may be thought of as follows. According 
to Section 2 (see the paragraph, preceding eq. (|2*7j) ). the volumes of the inner and outer 
spaces unite in sector IV. Consequently, the number of states is larger, than in sector II. 
The states in the negative-energy part of sector II are all occupied, whereas the newly 
added states are vacant. As a result, electrons in the Dirac see rearrange: they tend to 
leave the energy states with larger \e\ and keep to the border e = —m. The larger 
the more holes there are. 

The surface T_(e, R) is a maximum for every energy at R = —m' 1 . This fact 
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Figure 4. Transmission coefficient/positron distribution T_(e, R) (|51|l . i|56|) plotted 
in logarithmic scale against the negative energy ^ip- < for five negative values of 
the ratio Q, from left to right, Rm — — 1, — 2, — 3, — 0.2, — 0.1 The curves with 
\R\m < 1 (drawn solid) do not gather in the origin and show the threshold behaviour 
(E3 



is reflected in Figure 4: the curve for R = —m~ l occupies the extreme left position. 
The pair creation process runs most efficiently, when R = — m™ 1 , in other words, the 
distribution of produced particles is maximum at this value of R. The corresponding 
form of the distribution (jHTjl is 



(e, — m 



1-e 



27I-S. 

e p 



Not close to the threshold p — 0, this is especially simple: 



T_fe, 



-m 



e p . 



(62) 



(63) 



-i 



The same as in the positive-energy domain described above, the value \R\ = m 
discriminates two different situations. If \R\ > m _1 , R < 0, the border of sector IV 
is £ = —m, 7 > 0, and the bordering value of the distribution function (|51)l is zero: 

f_(-m,i?) = 0, \R\>m' 1 , (64) 

as is also seen in Figure 4: the curves with \R\ > m -1 gather in the origin T_ = 
0, e + m = 0. If \R\ < m _1 , R < 0, the border of sector IV is 7 = 0, e < —m. The 
threshold behaviour of the positron distribution is 

~ n(2) 3/2 VeR - 1 (l + coth^= 7r 

-m. (65) 



T_( £)J R) 
e<R- 1 < 



VI - m 2 i? 2 



This is zero in the threshold points e t hr = R 1 , which are the intersections between the 
curves with R > —mT l and the abscissa axis in Figure 4. By differentiating (J65|) over e 
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(66) 



near the threshold e = e t hr- The corresponding singularities are seen as peaks against 
the background of gentler and wider maxima in Figure 5. 
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Figure 5. Differential probability of positron creation - derivative of the transmission 
coefficient f|3T~fl . <|55|) dT_(e, i?)/d|e| plotted against the negative energy < 1 for 
five negative values of the ratio |HJ; from left to right, Rui — — 1, — 2, —3, —0.2, —0.1 
The curves with \R\m < 1 (drawn as solid) show the resonant threshold behaviour ll?)6l) 



We conclude that the absorption for < R < m _1 , e > m, and positron production 
for > R > — m _1 , e < — m are resonant near the thresholds of the corresponding 
processes. Note, that the condition \R\ < vnT 1 for the resonance to take place in the 
pair production differential probability coincides with the condition of the positivity of 
the norm discussed in Section 2. 

6. Discussion 

The existing theory of spontaneous e + e~-pair creation by a supercharged nucleus j^j, jH] 
does not solve the problem of singular interaction, but is satisfied with the pragmatically 
sufficient view, that the realistic nucleus has a finite size, which provides the regularizing 
cut-off of the potential near r = 0. The size of the nucleus core is taken to be about 
10~ 12 cm, which is an order of magnitude less than the only characteristic size in the 
problem, the electron Compton wave-length m~ l = 3.9 • 10 _11 cm. Consequently, all 
the results in that theory - including the value of the critical charge itself - are cut-off 
dependent and do not survive its removal. Besides, as distinct from ours, those are based 
on numerical calculations. These circumstances forbid a direct comparison, because in 
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our approach no cut-off is kept. Nevertheless we can point a certain correspondence 
between the spontaneous pair creation curves. 

The theory of [Sj, |H] points the critical value of the nucleus charge Z, the one for 
which the lowest electron level first sinks into the lower continuum e < —m and becomes 
a resonance, or a quasistationary state. The corresponding scattering amplitude acquires 
the Breit-Wigner shape. By considering a nonstationary problem, where the nucleus 
charge crosses its critical value in the course of time evolution, one establishes - with 
the use of the Fano theory [Hj - that the decay probability of the quasistationary state is 
determined by the same Breit-Wigner function. The quasistatinary state decays into a 
pair: the electron, which becomes a bound state localized near the centre, and the hole, 
which escapes to infinity and is interpreted as a free positron, subject to observation. 
(In our treatment, the created electron belongs to continuum of states, free near the 
center.) Thus, characteristic of the differential probabilities, obtained in this way, is 
a narrow Breit-Wigner resonance peak, its width depending on the cut-off. One may 
think, that the singular threshold behaviour seen in Figure 5 may be substituting for 
this resonance. 

As for the electron absorption by the supercritical nucleus, described in Section El 
and in [TO], this effect is unknown to the conventional theory. According to Figures 
1,2 the absorption runs very efficiently It might be observed if the heavy ion collision 
process is subjected to irradiation by a beam of electrons, or if the electrons produced 
in this process itself via, e.g., the two-photon mechanism are absorbed. 

There is another point, worth discussing. In our treatment above, the same as 
in the traditional treatment of the pair-production process by a supercharged nucleus 
[Sj, jH], the Dirac sea concept of the ground state and the hole theory of positrons 
were appealed to. It is known, however, that this concept cannot be done completely 
consistent, since the presence of a charge in the ground state remains not excluded 
within its scope. Nevertheless, this concept possesses a certain power of predictability 
and invokes useful analogies loaned from the solid state physics and theory of phase 
transitions. The interpretation of the produced positron distribution as resulting from 
the adding of a manifold of vacant states to the Dirac sea, proposed in Section |5] above, 
may be referred to as an example of this sort. 

In the meanwhile, the approach, developed in the present paper - see particularly 
Sectional where the singular problem is reduced to the barrier transmission/reflection 
on infinite axis - suggests a natural context for applying the theory of unstable vacuum, 
used in the literature to consider the Schwinger effect, i.e. particle production by an 
external electromagnetic field, pH], |19| . the Hawking radiation and Unruh's acceleration 
radiation jITj. This theory is based on second quantization and refers to nonequivalent 
Fock spaces, interrelated by the Bogoliubov transformation, the Bogoliubov coefficients 
being just the coefficients, that tangle asymptotes with positive and negative signs 
of momenta (in our case also of pseudomomenta) of the solutions to the equations 
of motion, prior to the second quantizaton. Leaving the detailed elaboration of the 
second-quantization programme of the singular problem to the forthcoming paper, we 
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now restrict ourselves to formulating a clue point, important for description of the 
electron-positron pair production. 

Let Oo Ut be operators, creating or annihilating states, which possess the wave 
function that behaves as one of the exponents (J7J) at the spatial infinity r — > oo. Besides 
the Fock space, spanned by the vectors created by repeatedly applying a^ nt to the 
out- vacuum |0 out ), <Vit|0out) = 0, define another Fock space, produced by the 
action on the true-vacuum |0 true ), a^ ue |0 true ) = 0, of the operator a^ ue , which 
creates states, whose wave function is given by the solution of the generalized eigenvalue 
problem, specified in Section El This eigenfuntion behaves as a linear combination of 
the two exponents ((Zj) at the spatial infinity r —>■ oo (also a linear combination of 
the two exponents (jHJ) or (JHUj) in the origin r — > or £ — > ±oo) - cf reference |14j . 
where analogous eigenproblem was explicitly solved for the Schrodinger equation with 
singular potential. To determine the mean number of particles in the true-vacuum 
state, suffices it to calculate the true-vacuum expectation value of the particle-number 
operator (O tr ue|aiut a out|0true)- This is expressed in terms of the coefficients in the linear 
combination of the exponents mentioned and corresponds to the count of particles by 
a remote observer. Simultaneously, another observer, who is placed in the origin - or is 
infinitely remote in the dilated-coordinate £-space, - would measure the mean number of 
particles (0t rue |ain~ a in|0true)) where are operators, creating or annihilating states 
with the wave function that behaves as one of the exponents (jHJ or (JHfij) in the origin 
r — > or £ — ► ±oo. The gases of particles, observed by the two observers, should be 
in mutual balance. The idea of balance is formally introduced into the theory, when 
we impose the non-Sturm boundary conditions ()15|) in the eigenvalue problem that 
interrelate the values of the wave function in the points r = and r = oo. 
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